Current-Driven Magnetic Excitations in Permalloy-Based Multilayer Nanopillars 
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We study current-driven magnetization switching in nanofabricated Ni84Fei6/Cu/Ni84Fei6 trilay- 
ers at 295 K and 4.2 K. The shape of the hysteretic switching diagram at low magnetic field changes 
with temperature. The reversible behavior at higher field involves two phenomena, a threshold cur- 
rent for magnetic excitations closely correlated with the switching current, and a peak in differential 
resistance characterized by telegraph noise, with average period that decreases exponentially with 
current and shifts with temperature. We interpret both static and dynamic results at 295 K and 
4.2 K in terms of thermal activation over a potential barrier, with a current dependent effective 
magnetic temperature. 

PACS numbers: 73.40.-c, 75.60.Jk, 75.70.Cn 



-a 
o 

> 
On 

m 
o 
m 
O 

-i— > 
c3 



i 

C 

O 

o 



X 



The effects of magnetic order on transport, e.g. gi- 
ant magnetoresistance, are well studied and widely used 
in technology. The discovery of predicted Q, Q current- 
induced magnetization precession 0, 0] and switching [j| 
stimulated an explosion of interest in the reverse effect 
of the current on magnetic order and magnetic dynam- 
ics. Besides fundamental interest in the physics of mag- 
netic systems driven far out of equilibrium, these phe- 
nomena hold promise for high-density memory appli- 
cations. With the number of theoretical papers grow- 
ing rapidly 0-0]' relatively few experimental facts 
are known, mostly for Co/Cu/Co multilayers, and, for 
nanofabricated samples, mostly at room temperature, 
T=295 K 

By studying Py/Cu/Py (Py=Permalloy=Nig4Fei6) 
nanofabricated trilayers (nanopillars), we are able for 
the first time to quantitatively compare current-driven 
switching at 295 K and 4.2 K. Much smaller crystalline 
anisotropy and magnetoelastic coefficients of Py elimi- 
nate the i rreg ular behavior at 4.2 K seen in Co/Cu/Co 
by others [19j and confirmed by us. Our data let us study 
temperature dependences, and establish general features 
of switching in both the hysteretic switching regime at 
low magnetic field H and the reversible switching regime 
at higher H. We emphasize two important findings. (1) 
In a new picture of the physics in the reversible regime, 
we isolate two different phenomena: (a) magnetic exci- 
tations occuring above a threshold current It, and ap- 
pearing either as a linear rise or peaks in differential re- 
sistance, dV/dl; and (b) the more well known reversible 
switching peak. We show that the latter arises from tele- 
graph noise switching of magnetization between paral- 
lel (P) and antiparallel (AP) orientations, and occurs 
when the dwell times in the P and AP states are ap- 
proximately equal, Tp ss tap- (2) Extending an idea of 
Ref. [2]]], we describe incoherent magnetic excitation by 
current in both hysteretic and nonhysteretic switching 
regimes in terms of an effective magnetic temperature 
T m , which can differ substantially from the lattice tem- 
perature T p h- Our model leads to a temperature depen- 



dence of the switching that differs from expectation for 
models 0] based on coherent current-driven excita- 
tion of the uniform precession [l| . 

Our samples were nanofabricated with a multistep pro- 
cess similar to that used by others . To minimize dipo- 
lar coupling between the Py layers, we used the geometry 
of Albert et al. {3| ! with the bottom (thicker) Py layer 
left extended, and the top Py layer and Cu spacer pat- 
terned into an elongated shape with typical dimensions 
130 x 60 nm. The bottom Py and the Cu spacer thick- 
nesses were always 20 nm and 10 nm, respectively (we 
give all thicknesses in nm). The patterned Py layer thick- 
ness was varied from 2 to 6 nm. We measured differential 
resistances, dV/dl, with four-probes and lock-in detec- 
tion, adding an ac current of amplitude 20 /xA at 8 kHz 
to the dc current I. Most of 12 Py/Cu/Py devices stud- 
ied had resistances R wl.5 fi, magnetoresistance (MR) 
«5% at 295 K, and MR«8% at 4.2 K. Positive current 
flows from the extended to the patterned Py layer. H is 
directed along the easy axis of the nanopillar. 

Figs. n^a,b) show the variations of dV/dl with I for a 
patterned Py(20)/Cu(10)/Py(6) trilayer at 295 K (la) 
and 4.2 K (lb), for H = 50 Oe (solid curves) and 
H = 500 Oe (dashed curves). The insets show the vari- 
ations with H for I — 0. At small H, the magnetization 
switches hysteretically to a higher resistance AP state 
at a large enough positive current lP~* AP = 7 S , and to 
a low resistance P state at negative lf p ~^ p . At larger 
H, the switching step turns into a nonhysteretic peak. 
Figs. [Tjc.d) show the switching diagrams at 295 K and 
4.2 K, extracted from data such as those in Figs.^a),(b), 
obtained both by varying I at fixed H and H at fixed I. 
Qualitatively, the 295 K data in Figs. GIa,c) are sim- 
ilar to those published previously for Co/Cu/Co flr^ |. 
As expected, both the reduced magnetization and ther- 
mal activation result in smaller switching currents and 
fields H S (I = 0) at 295 K. The slight ii-asymmetry in 
Figs.^c,d) is attributed to a combination of the current- 
induced Oersted field and sample shape asymmetry. 

The first important new feature of our data is the 
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FIG. 1: (a) Switching with current at 295 K in a 
Py(20)/Cu(10)/Py(6) trilayer (thicknesses are in nm). Solid 
line: H — 50 Oe, dashed line: H — —500 Oe. Arrows mark 
the scan direction. I t is the threshold current as defined in the 
text. Inset: MR dependence on H at I = 0. (b) Same as (a), 
at 4.2 K. (c),(d) Magnetization switching diagram, extracted 
from the current-switching at fixed H (solid symbols), and 
field-switching at fixed / (open symbols): (c) at 295 K, (d) at 
4.2 K. Downward triangles: AP to P switching, upward: P to 
AP switching. The reversible switching peaks are marked by 
coinciding upward and downward triangles. 
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FIG. 2: (a) Differential resistance of a 

Py(20)/Cu(10)/Py(3.5) trilayer for various H, at 4.2 K. 
Rp — 3.6 n, Rap = 4.0 Q, curves are offset for clarity. 
Dashed line follows the reversible switching peak. The onset 
of the peaks in Rp is marked by It, and 7 S is the P — » AP 
switching current at H = 0. (b) Relation between the 
P — > AP switching current I a at H = 0, T =4.2 K and the 
threshold current It (as defined in the text) for 12 samples 
with patterned Py layer thicknesses: 2 nm (squares), 3.5 nm 
(diamonds), 4.5 nm (triangles), 6 nm (circles). I 3 , It are 
multiplied by the resistances R of the samples. 



almost square shape of the hysteretic region at 4.2 K 
(Fig. ^d)). At 4.2 K, the switching currents do not 
change in the range —230 Oe< H < 230 Oe, beyond 
which the switching becomes reversible. Similarly, for 
—4 mA< I < 4 mA, H S (I) is independent of I, and at 
larger positive / the switching becomes reversible. 

The second important new feature is a threshold cur- 



rent It (labeled in Figs.^a,b) and Fig.|5fa)) seen at large 

H, with reversible switching. Most of our Py/Cu/Py 
and Co/Cu/Co samples showed a nearly linear rise of 
(dV/dl) p = Rp above It . Similar behavior was likely 
present in earlier Co/Cu/Co data at 295 K, but less ob- 
vious, and thus rarely noted [25|. In our smallest sample, 
(Fig. 0{a), dimensions estimated at 50 x 100 nm), the 
linear rise is resolved at 4.2 K into several peaks. They 
appear only to the left of the reversible switching peak, 
as the latter moves to higher I with increasing H (shown 
by a dashed line). We observed similar, less pronounced 
peak structures in some samples with smaller R. Re- 
markably, in all cases (e.g. Figs. da,b), Fig. I^a)), the 
onset of a linear rise or peak structures in Rp coincides 
with I S (H = 0). Within the uncertainty in determining 
the onset current It, we find I s — I t for all 12 samples 
studied (Fig. Gib)). This correlation between I s and It, 
combined with the abruptness of the changes with / in 
Fig.^d), point to the importance of It- Moreover, upon 
following the data up to H — 4 kOe (not shown) , we find 
It to be almost independent of H (increasing by only 
10-40% for different samples). 

We identify It as a threshold for magnetic excitation 
of the patterned Py layer. An excitation threshold was 
previously seen as a peak in the differential resistance of 
point-contacts on extended magnetic multilayers 0, 0] . 
In nanopillars, this threshold behavior is modified into 
a linear rise or peak structures. The reversible switch- 
ing peak appears at higher H , and evolves faster with H 
(dashed line in Fig.[3{a)) than the peaks we attribute to 
magnetic excitations. Time-resolved measurements of re- 
sistance, at I and H close to the reversible switching peak 
(Fig. 01a)), show that the reversible switching is charac- 
terized by telegraph noise with random distribution of 
Tp, tap- Similar slow telegraph noise was reported in 
a point contact at one current and at the transition 
point from hysteretic to reversible switching in Co/Cu 
nanopillars at 295 K, at I or H fixed [Tij. 

The data of Fig. EJb,c), and the finding that the re- 
versible switching peak is a consequence of telegraph 
noise and occurs at rp w tap, represent our third new 
experimental result. Fig. E^b) shows that, when both I 
and H are increased so as to hold rp = tap, the aver- 
age period of the telegraph noise decreases exponentially 
with similar slopes at 295 K and 4.2 K, down to the 
1 MHz bandwidth limit of our setup. Fig. UJc) shows 
that the variations of tp,tap with I have similar forms 
at 295 K and 4.2 K. These data (as well as similar data 
for Co/Cu at 295 K [3 Hi]) show that the reversible 
switching peak does not indicate abrupt onset of a new 
physical process. At small /, rp >> tap, and the av- 
erage sample resistance is close to Rp, while at large 

I, Tp « tap, and the average resistance is close to 
Rap- The reversible switching peak appears in dV/dl 
at Tp tap, due to the exponential variation of rp, tap 
with /. Thus, it merely reflects the current-dependent 
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FIG. 3: (a) Time traces of the Py(20)/Cu(10)/Py(6) sam- 
ple static resistance R = V/I at 7 = 4.4 mA, T =295 K 
at different H. (b) Current dependence of the average tele- 
graph noise period for Py(20)/Cu(10)/Py(6) (squares) and 
Py(20)/Cu(10)/Py(3.5) (circles) samples. Open symbols: 
T =4.2 K, filled symbols: T =295 K. H was adjusted so that 
the average dwell times in AP and P states were equal, (c) 
The dependence of the P (downward triangles) and AP (up- 
ward triangles) dwell times on 7, for a Py(20)/Cu(10)/Py(6) 
sample. Open symbols: T =4.2 K, H = -335 Oe, filled 
symbols: T =295 K, 77 = -120 Oe. (d) Dependence of P 
(downward triangles) and AP (upward triangles) dwell times 
on H at I = 4.4 mA, T = 295 K. 



telegraph noise statistics; I t is the only fundamental cur- 
rent related to magnetic excitations. The amplitude and 
inverse width of the reversible switching peak are pro- 
portional to d[(rAp — Tp) I '(tap +Tp)]/dI. In Fig. Eta), 
they are correlated with the positions of the peaks that 
we attribute to magnetic excitations. For example, at 
200 Oe, the reversible peak is on top of the magnetic 
excitation peaks, and is significantly taller and narrower 
than at higher 77. More detailed data and analysis will 
be presented elsewhere |2fj|. 

We interpret both our static and dynamic results con- 
sistently in terms of a model involving thermally acti- 
vated switching over a one-dimensional potential bar- 
rier, with an effective current-dependent temperature 
T m (7) 21], and a threshold current, I t . T m can dif- 
fer from the lattice temperature T p h in confined ge- 
ometry, where the magnetic energy relaxes significantly 
more slowly than the highly excited individual magnetic 
modes [2^ |. Formally, a potential energy barrier is in- 
consistent with the current-driven torque in the Landau- 
Lifshitz- Gilbert equation Li and Zhang have 
shown that the work performed by the torque leads to 



a current-dependent effective switching barrier, but with 
T m = Tph, which we shall see is inconsistent with our 
data of Fig.Etb,c). 

Our observation that 7 t varies only weakly with 77 up 
to 4 kOe provides information about whether its ori- 
gin is quantum mechanical Q or classical 0. In the 
quantum model, 7 t follows from A/i = huj, where A/z 
is the current dependent difference in chemical poten- 
tials of spin-up and spin-down electrons due to spin ac- 
cumulation, and lu is the magnpn frequency. At small 
if, It oc oj oc \/l + H/ H a |22j) where H a is the in- 
plane anisotropy field. In our samples, H a pa 300 Oe, 
giving much too rapid an increase of It with 77. In 
the classical treatment based on the Landau-Lifshitz- 
Gilbert equation, I t is set by the balance between the 
torque induced by the spin polarized current and Gilbert 
damping. The H dependence of It is much weaker 
7 t (77) w J t (0)[l + H/(2itM)] [|, where M is the mag- 
netization, M k 880 Oe for Py. This model predicts 
i!(a)dH ~ 0-2 (kOe) -1 , not too far from the measured 
values of 0.03-0.1 for different samples. 

To analyze our data quantitatively, we determine the 
dwell times tp.ap by 



tp.ap 



1 
9. 



exp 



U 



P,AP 



kT, 



P,AP 



and approximate T m by the heuristic relation 
K(I-It) for I > I t . 



ph 



(1) 



(2) 



Here K is a constant, 9 pa 10 7 s 1 |29j is the effective at- 
tempt frequency, Up^ap is the potential barrier height for 
switching from the P or AP state, and T£ P (I) T%(I) 
when 7^0. The switching barriers U p,ap depend on 
7 only through the variation of the magnetization with 
temperature T m . In the P state, at I > It, magnetic exci- 
tation leads to increase of T^, as illustrated in Fig.QJa). 
A thermally activated transition into the AP state occurs 



at kT^ 
time t 



Up 



U^-, based on the data acquisition 



/n(t exp n) ~ 16 

exp of 1 second per point. The conditions for mag- 
netic excitation are not satisfied in the AP state, so the 
magnetic system cools to T™^ w T p h and at 77 < H s 
becomes trapped in this state. At 77 > H s (Fig. 0Jb)), 
kTph > ■=^ £ , i.e. the AP — > P transition is also ther- 
mally activated, leading to telegraph noise at 7 > 7 S , 
77>77 s . 

The 4.2 K switching diagram of Fig. ^d)) is consistent 
with this model. It reflects the weak variations of It 
with 77, and of T m (and thus 77 s ) with 7 below 7 t . We 
attribute the rounding of the 295 K diagram at 7 < to 
enhancement of thermal fluctuations of magnetization by 
current, not included in the heuristic Eq. J5J. Increasing 
H decreases Uap , so that even these weak excitations can 
activate the AP^P transition. This process is described 
in different terms in Ref. |l5j. 
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FIG. 4: (a) Schematic of current-driven hysteretic magneti- 
zation switching at H — 0. Dashed lines indicate T m . (b) 
Schematic of telegraph noise at H > H s . 



The dependencies of nonhysteretic switching on /, H, 
and T shown in Fig. Etb-d) agree with Eqs. QJ and J2J. 
As / increases with H fixed, rp decreases exponentially 
because of the increase of with I > I t . An increase of 
tap with / may be evidence for cooling of the magnetic 
system of the patterned layer in the AP state at / > 0, 
through a magnon absorption mechanism inverse to the 
magnon emission in the P state. At T p h,=4.2 K, the re- 
quirement for telegraph noise (Uap < 16fcT p ^) is satisfied 
with much smaller Uap than at 295 K. Hence a slight de- 
crease in T^ P with I has a stronger effect on tap at 4.2 K 
than at 295 K, as seen in Fig. Etc). The dependencies in 
Fig. Eld) follow from dlnT a p ZJ AP = ,J, A ^ dU ^ p . Since 



dH 



kT„ 



dH 



> T£ , Tp varies with H more slowly than tap- 
Based on the data of Fig. [31 we estimate the dependence 
T m {I) numerically. Some of our Py(20)/Cu(10)/Py(3.5) 
and all of the Py(20)/Cu(10)/Py(2) samples were su- 
perparamagnetic at 295 K. Assuming linear variation of 
the switching barrier with nanopillar thickness, we esti- 
mate the switching barrier Up } ap(H = 0) » 0.7 eV for 
Py(20)/Cu(10)/Py(6) samples. We use an approximate 
dependence U(T m ) = U ^l -T m /T c , where T c = 800 K 
is the Curie temperature for Py. Using d\u(Tp)/dI 
1 1 (mA)" 1 from Fig. Etc), we obtain from Eq. 

K = S£ « 400 K/mA. Both the 295 K (filled sym- 
bols in Fig. Etc)) and 4.2 K data (open symbols) are ap- 
proximately consistent with this estimate. In contrast, 
starting from a current-dependent effective barrier with 
T m = T ph pi 01, one would predict dln( g (J)) « jjjfi-, 
where I c is the switching current at T = 0. This strong 
dependence on T p h disagrees with the similarity of the 
4.2 K and 295 K data. 

To summarize, our major new experimental results on 
Py/Cu/Py nanopillars are: i) a square switching dia- 
gram at 4.2 K; ii) an onset current It (closely related to 
the hysteretic switching current I s ) for a linear rise of 
dV/dl in larger samples or a series of peaks in smaller 
ones; iii) reversible switching, characterized by telegraph 
noise with rate both increasing exponentially with / and 
shifting with temperature. The reversible switching peak 
in dV/dl occurs when the dwell times in the P and AP 
states are approximately equal. We are able to explain 



the hysteretic switching behavior at both 295 K and 
4.2K, and variations of telegraph noise with /, H , and 
T, by means of a threshold current for the onset of mag- 
netic excitations and thermally activated switching with 
a current-dependent magnetic temperature. 
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